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The nonlinear coefficients defining the mean electromotive force (i.e., the nonlinear turbulent 
magnetic diffusion, the nonlinear effective velocity, the nonlinear K-tensor, etc.) are calculated 
for an anisotropic turbulence. A particular case of an anisotropic background turbulence (i.e., 
the turbulence with zero mean magnetic field) with one preferential direction is considered, ft is 
shown that the toroidal and poloidal magnetic fields have different nonlinear turbulent magnetic 
diffusion coefficients. It is demonstrated that even for a homogeneous turbulence there is a nonlinear 
effective velocity which exhibits diamagnetic or paramagnetic properties depending on anisotropy 
of turbulence and level of magnetic fluctuations in the background turbulence. The diamagnetic 
velocity results in the field is pushed out from the regions with stronger mean magnetic field, 
while the paramagnetic velocity causes the magnetic field tends to be concentrated in the regions 
with stronger field. Analysis shows that an anisotropy of turbulence strongly affects the nonlinear 
turbulent magnetic diffusion, the nonlinear effective velocity and the nonlinear a-effect. Two types 
of nonlinearities (algebraic and dynamic) are also discussed. The algebraic nonlinearity implies a 
nonlinear dependence of the mean electromotive force on the mean magnetic field. The dynamic 
nonlinearity is determined by a differential equation for the magnetic part of the a-effect. It is shown 
that for the att axisymmetric dynamo the algebraic nonlinearity alone (which includes the nonlinear 
a-effect, the nonlinear turbulent magnetic diffusion, the nonlinear effective velocity, etc.) cannot 
saturate the dynamo generated mean magnetic field while the combined effect of the algebraic and 
dynamic nonlinearities limits the mean magnetic field growth. 

PACS numbers: 47.65.+a; 47.27.-i 



I. INTRODUCTION 



equation 



dB/dt = Vx(VxB + 5-7 ? VxB), 



(1) 



Generation of magnetic fields by turbulent flow of 
conducting fluid is a fundamental problem which has 
a large number of applications in solar physics and as- 
trophysics, geophysics and planetary physics, etc. In 
recent time the problem of nonlinear mean-field mag- 
netic dynamo is a subject of active discussions (see, e.g., 
§ | | | | | |, | | 0). It was suggested in Q 
that the quenching of the nonlinear a-effect is very strong 
and causes a very weak saturated mean magnetic field. 
However, the later suggestion is in disagreement with ob- 
servations of galactic and solar magnetic fields (see, e.g., 
H H PJR16L01) and with numerical simulations 
(see, e.ffTiMllP 

Saturation of the dynamo generated mean magnetic 
field is caused by the nonlinear effects, i.e., by the back 
reaction of the mean magnetic field on the a-effect, tur- 
bulent magnetic diffusion, differential rotation, etc. The 
evolution of the mean magnetic field B is determined by 
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where V is a mean velocity (e.g., the differential ro- 
tation), rj is the magnetic diffusion due to the electri- 
cal conductivity of fluid. The mean electromotive force 
£ = (u x b) in an anisotropic turbulence is given by 

£, = a l3 B J + (V cff xB) l -r !l3 (VxB) 3 

-K ijk (dB) jk -[Sx(VxB)]i (2) 

(see |H H), where (dB\ 3 = (1/2)(V 1 B J + Vj-Bj), u 
and b are fluctuations of the velocity and magnetic field, 
respectively, angular brackets denote averaging over an 
ensemble of turbulent fluctuations, the tensors ay and 
rjij describe the a-effect and turbulent magnetic diffu- 
sion, respectively, V cff is the effective diamagnetic (or 
paramagnetic) velocity, Kijk and 8 describe a nontriv- 
ial behavior of the mean magnetic field in an anisotropic 
turbulence. Nonlinearities in the mean-field dynamo im- 
ply dependencies of the coefficients (a^ ;, rjij, V off , etc.) 
defining the mean electromotive force on the mean mag- 
netic field. The a-effect and the differential rotation are 
the sources of the generation of the mean magnetic field, 
while the turbulent magnetic diffusion and the K-effect 
(which is determined by the tensor Kijk) contribute to 
the dissipation of the mean magnetic field. 
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In spite the nonlinear a-effect was under active study 
(see, e.g., & R]), the nonlinear turbulent magnetic dif- 
fusion, the nonlinear K-effect, the nonlinear diamagnetic 
and paramagnetic effects, etc. are poorly understood. 

In the present paper we derived equations for the 
nonlinear turbulent magnetic diffusion, the nonlinear 
effective velocity, the nonlinear K-effect, etc. for an 
anisotropic turbulence. The obtained results for the 
nonlinear mean electromotive force are specified for an 
anisotropic background turbulence with one preferen- 
tial direction. The background turbulence is the turbu- 
lence with zero mean magnetic field. We demonstrated 
that toroidal and poloidal magnetic fields have different 
nonlinear turbulent magnetic diffusion coefficients. It is 
shown that even for a homogeneous turbulence there is 
a nonlinear effective velocity which can be a diamagnetic 
or paramagnetic velocity depending on anisotropy of tur- 
bulence and level of magnetic fluctuations in the back- 
ground turbulence. 



II. THE GOVERNING EQUATIONS 

In order to derive equations for the nonlinear turbu- 
lent magnetic diffusion and other nonlinear coefficients 
defining the mean electromotive force we will use a mean 
field approach in which the magnetic, H, and velocity, 
v, fields are divided into the mean and fluctuating parts: 
H = B + b, v = V + u, where the fluctuating parts 
have zero mean values, V = (v) = const, and B = (H). 
The momentum equation and the induction equation for 
the turbulent fields u and b in a frame moving with a 
local velocity of the large-scale flows V are given by 



db 
~dt 



VP' 



f 



-[b x (V x B) + B x (V x b)] 

P P-oP 

-T + Mu + F/p , 



V x (u x B - 77V x b) + G 



(3) 
(4) 



and Vu = 0, where P' are the fluctuations of the hy- 
drodynamic pressure, F is a random external stirring 
force, v is the kinematic viscosity, 77 is the magnetic dif- 
fusion due to the electrical conductivity of fluid, p is the 
density of fluid, fiQ is the magnetic permeability of the 
fluid, the nonlinear terms T and G are given by T = 
((u- V)u) - (u- V)u+ [<b x ( V x b)) - b x (V x b)]/( M o/o), 
and G = Vx(uxb — (ux b)). We consider the case 
of large hydrodynamic (Re = I^uq/v 3> 1) and magnetic 
(Rm = IqUq/t] S> f) Reynolds numbers, where uq is the 
characteristic velocity in the maximum scale of turbu- 
lent motions. 



The procedure of the derivation of equation for 
the nonlinear mean electromotive force 



The procedure of the derivation of equation for the 
nonlinear mean electromotive force is as follows (for de- 



tails, see Appendix A). 

(a) . By means of Eqs. (||) and (^) we derive equations 
for the second moments: 

/«(k,R) = y< Uj (k + K/2) Uj (-k + K/2)) 

x exp(iK-R) dK = /ji(-k, R) , (5) 

/iy(k,R) = /"{fciOc + K/2)6j-(-k + K/2)) 

x exp(iK-R) (IK/ nop = hji(—k, R) , (6) 

S«(k,R) = /<&i(k + K/2K-(-k + K/2)) 

x exp(iK-R) dK . (7) 

where R and K correspond to the large scales, and r 
and k to the small ones, i.e., R = (x + y)/2, r = 
x-y, K = k!+k 2 , k=(ki-k 2 )/2. 

(b) . We split all correlation functions (i.e., fy, hij, gy) 

into two parts, e.g., hij 



AN) 



where the 



tensor — [/iy(k,R) + hy(— k,R)]/2 describes the 



nonhclical part of the tensor and h 



[hij(k, R) 



hij(— k, R)]/2 determines the helical part of the tensor. 
Such splitting is caused, e.g., by different times of evolu- 
tion of the helical and nonhclical parts of the magnetic 
tensor. In particular, the characteristic time of evolution 
of the tensor is of the order tq = Iq/uq, while the 
relaxation time of the tensor h\ - ' is of the order of ToRm 
(see, e.g., (l§ || || ||). 

(c). Equations for the second moments contain higher 
moments and a problem of closing the equations for the 
higher moments arises. Various approximate methods 
have been proposed for th e solution of problems of this 
type (see, e.g., 
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|2q]). The simplest procedure is 
the r-approximation, which is widely used in the the- 
ory of kinetic equations. For magnetohydrodynamic tur- 
bulence this approximation was used in ||29f (see also 
0, |3l]]). In the simplest variant, it allows us to express 
the third moments in terms of the second moments: 



Mu - M. 



(0) 



R 



13 

(AO 
ij 



-(/«-/i3V(*o 



hS N) )/r(k) , 



(8) 
(9) 
(10) 

eqi 



where My, Rij, Cij are the third moments 
tions for fij,hij and gij, respectively (see Eqs. ( A3)- 
( |A"e| ) in Appendix A). The superscript (0) corresponds 
to the background magnetohydrodynamic turbulence (it 
is a turbulence with zero mean magnetic field, B = 0), 

h^ N ^ is the nonhelical part of the tensor of magnetic fluc- 
tuations of the background turbulence, and r(fc) is the 
characteristic relaxation time of the statistical moments. 
We applied the r-approximation only for the nonhelical 

part hj*p of the tensor of magnetic fluctuations because 

(s) 

the corresponding helical part h\- is determined by an 
evolutionary equation (see, e.g., B @, Il5l p3|, p3L EM |33] 
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and Section III-C). We took into account here magnetic 
fluctuations which can be generated by a stretch-twist- 
fold mechanism when a mean magnetic field is zero (see, 
e.g., |3^, Q). This implies that nf? ^ 0. In inertia 
range of background turbulence i?y (B = 0) — and 
dj (B = 0) = 0. We also took into account that the cross- 
helicity tensor gij for B = is zero, i.e., <7ij(B = 0) = 0. 

The r-approximation is in general similar to Eddy 
Damped Quasi Normal Markovian (EDQNM) approx- 
imation. However some principle difference exists be- 
tween these two approaches (see J2(| |28|]). The EDQNM 
closures do not relax to equilibrium, and this procedure 
does not describe properly the motions in the equilib- 
rium state in contrast to the r-approximation. Within 
the EDQNM theory, there is no dynamically deter- 
mined relaxation time, and no slightly perturbed steady 
state can be approached p6| . In the r-approximation, 
the relaxation time for small departures from equilib- 
rium is determined by the random motions in the equi- 
librium state, but not by the departure from equilib- 
rium p6| . We use the r-approximation, but not the 
EDQNM approximation because we consider a case with 
lo\V B 2 \/ [i <C (pu 2 ). As follows from the analysis by [2f| 
the r-approximation describes the relaxation to equilib- 
rium state (the background turbulence) much more ac- 
curately than the EDQNM approach. 

In this study we consider an intermediate nonlinearity 
which implies that the mean magnetic field is not enough 
strong in order to affect the correlation time of turbulent 
velocity field. The theory for a very strong mean mag- 
netic field can be corrected after taking into account a 
dependence of the correlation time of the turbulent ve- 
locity field on the mean magnetic field. 

(d). We assume that the characteristic time of varia- 
tion of the mean magnetic field B is substantially larger 
than the correlation time r(fc) for all turbulence scales. 
This allows us to get a stationary solution for the equa- 
tions for the second mome n ts fa , hij and gij . Using these 
equations [see Eqs. ( Al3| )-( A20) in Appendix A] we cal- 
culate the electromotive force £j(r = 0) = J Si (k) dk, 

where £ t (k) = (l/2)e. 
result is given by 



i(«$?(k,R) 



5 S(-k,R)).The 



£i(r = 0) = aijBj + /..,;, //,> , 
where Bij — dBi/dRj, 



(11) 



an = * / r(l + ^)- 1 £ „„„fc J (/(° m s )-/ 1 ( l f))dk,(12) 



.Oc. h ./) Wi die 



(13) 



(for details, see Appendix A), fcjj = kikj/k 2 , /il™ = 

h>nm ) +ip(l + 2ip)- 1 (fi°m ) -h < nm ,> ), s ijk is the Levi-Civita 
tensor, and ip = [(/3-k)u r/2] 2 , fa = AB l / (u Qy /2jI^p), 

fij N ^ and flj^ describe the nonhelical and helical ten- 
sors of the background turbulence. 



(e). Following to Q] we use an identity Bj, k = 
(dB)jk — £jfcz(VxB);/2 which allows us to rewrite Eq. 
(lilf) for the electromotive force in the form 



Si = aijBj + (UxB), - j7y(VxB)j - n ijk (dB) jk ,(14) 
where 

Oij(B) = (dij +ciji)/2 , C4(B) = EkjiCHj/2 ,(15) 
Kij k (B) = -(hjk + b ik j)/2 . (16) 

B. The model for the background turbulence 

For the integration in k-space in Eqs. (JlJ) and jl^ ) 
we have to specify a model for the background turbulence 
(i.e., turbulence with zero mean magnetic field). We as- 
sume that the background turbulence is anisotropic and 
incompressible. The second moments for turbulent ve- 
locity and magnetic fields of the background turbulence 
are given by 

r Cii (k) = (5/4){P y -(fc)[(2/5)77 T l) (k)- Ai W(k)fc nm ] 
+2[% M W (k)k nm + 4 a) (k) - M&!(k)*W 



-*w45(k)]} 

(see 0), where = f^ N ^ when a = v, and c 



(17) 

(ON) 



when a = h, and ry^ (k) = r/pp JV; (k), fj^'Qc) = 
Th^p p N ^ (k), Pij (fc) = Sij — kij, S mn is the Kronecker tensor. 
The anisotropic part of this tensor /Xmn (k) has the prop- 
erties: jUmn(k) = /inm(k) and /Xp^(k) = 0. Inhomogene- 
ity of the background turbulence is assumed to be weak, 
i.e., in Eq. ( p7[ ) we dropped terms ~ 0[V (rj^ ; My )L 

where = tqUq/3, 77^ = ro6§/3/ioP and &o is the char- 
acteristic value of the magnetic fluctuations in the back- 
ground turbulence. To integrate over k in Eqs. ( |l2| ) and 
( |l3| ) we use the Kolmogorov spectrum of the background 

turbulence, i.e., Tfpp N \\i) = rff\p(k), Th p p N '' '(k) = 
rf^ip(k) and jti£m(k) = ^4n„(R)^?(fc)/3, where tp(k) = 
(■Kk 2 k )-\k/k )- 7 / 3 , r{k) = 2r (fc/fc )- 2 / 3 , k = Iq 1 . 
We take into account that the inertial range of the turbu- 
lence exists in the scales: Id < r < Iq. Here the maximum 
scale of the turbulence Iq <C Lb, and Id = lo/Re 3 ^ 4 is 
the viscous scale of turbulence, and Lb is the character- 
istic scale of variations of the nonuniform mean magnetic 
field. 

In the next section we present results for the nonlinear 
coefficients defining the mean electromotive force. 



d h )t 



III. NONLINEAR COEFFICIENTS DEFINING 
THE MEAN ELECTROMOTIVE FORCE 

The procedure described in Section II (see also, for 
details Appendix A) allows us to calculate the nonlin- 
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ear turbulent magnetic diffusion tensor, the nonlinear k- 
tcnsor, the nonlinear a-tensor and the nonlinear effective 
drift velocity. 



B. 



The hydrodynamic part of the nonlinear 
a-tensor 



Using Eqs. (|f 2j) and (M) we get 



A. Nonlinear turbulent magnetic diffusion tensor 
and nonlinear K-tensor 



The general form of the turbulent magnetic diffusion 
tensor %-(B) contains all possible tensors: #y, /iy fyj 



and their symmetric combination /n 



(«) 



[see Eq. ( A50 ) in Appendix A], where = (3 % f3j/Pi 2 , 
fii = &BiJ (uq\/2iiqp). For an isotropic background tur- 
bulence (when fj,\j' = 0) the turbulent magnetic diffusion 
tensor r?y(B) is given by 

% (B) = SijiAiiy/lp)^ + [Ai(0) - Ax(V2(3)]r}^} 
+ 1 (3 ij A 2 {0){r 1 ^ +r,P), (18) 



where the functions Ak(P) are defined in Appendix B. 
For P < 1 Eq. @ reads 



% -(B) = %[ 7? W-(2/3 2 /5)(2^ ) -^ ) )] 



-(2/5)AA-(^ ) +r?? ) ) 



(19) 



and for /3 ^> 1 it is given by 

% -(B) = (Stf/IO/?)^^-^ 
-/3y(^+4 h) )}. 



W (^- 



1)] 



(20) 



The mean magnetic field causes an anisotropy of the tur- 
bulent magnetic diffusion tensor which is determined by 
the tensor Magnetic fluctuations of the background 
turbulence contribute to the turbulent magnetic diffusion 
tensor 77^ (B) in the nonlinear case, ft follows from Eq. 
( p0| ) that for [3 ^> 1 the tensor rjy oc 1//3. 

The K-tensor describes a nontrivial behavior of the 
mean magnetic field in an anisotropic turbulence. For an 
isotropic background turbulence the K-tensor vanishes in 
spite of an anisotropy caused by the mean magnetic field. 
For an anisotropic background turbu lence a general form 
of the K-tensor is given by Eq. ( A51 ) in Appendix A. For 
P -C 1 this tensor is given by 



f(v) f(h)} 
J ijk ' ^ijk 



-(l/6)(3L^ + L^) + (1/7)^(5^ 



(21) 



and for /3 3> 1 it reads 



where L 
.(«)/ 



(a) 



= -(7r/16/3)(V2-l)[L (,;) 



ijfc 



f 00 



(a) 



ijk l\ ' 

(a) 

£ikn^ n j 



(22) 



and 



AT. 



Cr( a ) 



Hnp( £ ijnP P k + £iknP P j). Note that for /3 > 1 the ten- 
sor oc 1//3. The K-tensor contributes to the turbu- 
lent magnetic diffusion of the toroidal and poloidal mean 
magnetic fields (see Section V). 



<(B,R) = 



l + V(B,k) 



(23) 



where hereafter aff (0, k, R) = a^(B = 0,k,R). Anal- 
ysis in j?]. H^] shows that a form of the tensor oiy (0, k, R) 
in an anisotropic turbulence can be constructed using the 



tensors fey, fcy TOn and 14 



y , where 



is the anisotropic part of the hydrodynamic contribution 
of the a-tensor. Thus we use the following model for the 
tensor a^(0, k, R) 



ag>(0,k,R) = {2a^'(R)k ij + 5ek ijmn v mn {R) 

+ (1 - e)[v ip (R)k pj + u jp (K)k pi }}ip(k)/2 , (24) 

where the parameter e describes an anisotropy of the he- 
lical component of turbulence and it changes in the in- 



» 



terval: < e < 1. Here a£°(0,R) = J a^ ; (0,k,R)dk 



(v). 



(R)5i. 



(R), and 4 U) (R) = (l/3)4p ; (0,R), the 



anisotropic part Vij (R) of the hydrodynamic contribution 



of the a-tensor has the properties: fj. 
Substituting 



Vji and 



0. 

Eq. (|24j) into Eq. (|23j), and using identities 
30) we obtain the nonlinear dependence of 



(|A2S| ) and 

the hydrodynamic part of the a-effect on mean magnetic 
field: 



(B) = (l/2){^[2(A 1 (/3)+A 2 (/3))a^ 

+ (1 - e)A 2 {fi)u p + 5e(C 2 (/3) + 3C 3 (/3))^] 

+^[(l-e)(2A 1 (/3) + J 4 2 (/3)) 

+ 10e(Ci(/3)+C 2 (/3))]} , (25) 



where ^(R) = v mn (R)l3 nm and the functions Ck(/3) are 
defined in Appendix B. For e = Eq. ( |BJ ) coincides with 

that derived in . The asymptotic formulas for a-^ for 
/? -C 1 and (3 > 1 are given by Eqs. flA53j ) and (|A56|) in 
Appendix A. 



C. The mean electromotive force and the nonlinear 
magnetic a-tensor 



Using Eqs. ( |A40[ ), (JA43J) , ( |A46| ) we calculate the elec- 
tromotive force E 



-K ijk (dB) jk , (26) 

where the nonlinear effective drift velocity V cff = U + 
V (JV) , and the velocity J7j(B) = -(l/2)e imn o m „ = 
-(l/2)VpA^f (see §), the velocity is given 

by Eq. ( A44| ) , the tensor of turbulent magnetic diffusion 
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r/ij is given by Eq. (A50), the tensor Kyj, is determined 



anisotropic part \x -J ^ of the tensor ^J" 1 (B = 0) is given by 



(M) is defined in Eqs. ( |A25| ) 



by E q. (|A5l|) , the tensor A y 
and (A33). In the kinematic dynamo the effective drift 
velocity (turbulent diamagnetic velocity) is caused by an 
inhomogeneity of turbulence. The effective drift veloc- 
ity U(B) is determined by the tensor ay and is due to 
an induced inhomogeneity of turbulence caused by the 
nonuniform mean magnetic field. This implies that the 
nonuniform mean magnetic field modifies turbulent ve- 
locity field and creates the inhomogeneity of turbulence. 
The effective velocity V^^B) is determined by tensor 
bijk and is caused by the nonuniform mean magnetic field. 
The a-tensor is determined by the hydrodynamic and 



0) 
My 
» 



(ey-(l/3)<5y), and/4? = l4„ Pnj + Pinl*. 



(v) 



[(eifo + e j /3 l )(e-$) - (2/3)/%], where ej^ is a de- 



gree of of an anisotropy of the turbulence, and /A 



(v) 



(l/2)j2 ( pj = e^'[(e-0) 2 -1/3]. It follows from these equa- 



tions that ^(B = 0) = - (1/3)4^] + eye 



(v) 



.(«) 



Now we take into account that the components r/xx (B = 

0), rjy V J (B = 0) and rfzz (B = 0) are positive. This yields 

—3/2 < /rf>p < 3. The equations for the correspond- 
ing magnetic tensors are obtained from these equations 
after the change v — > h. For the magnetic fluctuations we 



magnetic contributions, i.e., ay (B) 
with 



» 



(B)- 



(B) also obtain that -3/2 < ef'/rjf' < 3 



a^(B) = 4 h) (B)$(/3),yy 



(27) 



(see |tJ ) , where the tensor af^ (B) is determined by Eq. 
@, the function $(/3) = (3//3 2 )[l - arctan(/3)//3], and 

the magnetic part ag^(B) of the a-effect is determined 
by the dynamic equation 



da, 



(h) 



dt 



(h) 

T 



'flux) 

£(B)B 



For galaxies, e.g., the preferential direction e is along 
rotation (which is parallel to the effective gravity field). 
For the axisymmetric afi-dynamo and large magnetic 
Reynolds numbers the toroidal magnetic field is much 
larger than the poloidal field. Therefore, the value e-/3 
is very small and can be neglected because (3 is approxi- 
mately directed along the toroidal magnetic field. 

Thus, the nonlinear coefficients defining the mean elec- 
tromotive force in a turbulence with one preferential di- 
rection are given by 



r?y (B) = MrjSij + M e e i:j + Mpfiij , (30) 



(28) 



V cff (B) = B 



~ 2 [M^VB 2 



M v 2) e(e-V)B 2 



(31) 

K ijk (B){dB) jk = M K [ex(e-V)B] i (32) 
(see |§§§§§), where W t = Clj V 3 is the 

velocity which depends on the mean fluid velocity V (see Appendix C), where we assumed that e-/3 = 0, the 
(for an isotropic turbulence the tensor cy = 5y and 
for an anisotropic turbulence with one preferential direc- 
tion, say in the direction e, the tensor cy = (23/30)(5y + 



functions M n , M e , Mp, M K , My and My' are given 
by Eqs. (0)-(jC|) in Appendix C. The tensor r?y(B) 



r(2) 



(7/10)ejej, see 25 ); the flux 



,o»(B)M jfm* 

P V?7^ ) (B = 0)/xoP 



F flu 



(29) 



contains three tensors <5y, ey and /3y since here there 
two preferred directions, along the vectors e and B. 

Now we consider the hydrodynamic part of the a- 
effect for an anisotropic background turbulence with 



is related with the flux of the magnetic helicity and is in- 
dependent of the mean fluid velocity V || (see also [[l0| ) , 
and T ~ roRm is the characteristic time of relaxation of 
magnetic helicity. The asymptotic formulas for for 



one preferential direction. The tensor a^(B 
this case can be rewritten in the form a)" (B 



r O) 



j3 < 1 and > 1 are given by Eqs. (A54) and (A57) in 
Appendix A. 



IV. ANISOTROPIC BACKGROUND 
TURBULENCE WITH ONE PREFERENTIAL 
DIRECTION 

Now we consider an anisotropic background turbulence 
with one preferential direction, say along an unit vec- 
tor e. Thus the tensor iffi (B = 0) 



0) in 

0) = 

a o&ij + v %j = [ a o"^ - (l/ 3 ) £ a]5y + £ Q ey, where e a is 
a degree of an anisotropy of the a-tensor. Thus, the 
anisotropic part ty is given by z/y = £ Q (ey — (l/3)<5y). 
The electromotive force contains the tensor ay in the 
form ctijBj. Thus, fy/3j = — (l/3)e Q [/?i — 3(e-/3)e,] and 
up = -(l/3)e a [l - 3(e-/3) 2 ]. Using Eq. (§|) we obtain 
the hydrodynamic part of the a-tensor in an anisotropic 
background turbulence with one preferential direction 



4 } (B) 



^•{[^iCS) + A 2 (p)][a^ - (l/3)e Q (l - e)] 
-(5/6)e a e[2Ci(/9) + 3C 2 (/3) + 3C 3 (fi)]} 
= $ a (P)6ij, (33) 



by % (B - 0) - rj T dij + /x y - % dy 



TViVj) is given 

( , 

£ u e? 



where we assumed that e-/3 = 0. For e / the tensor 



» 



(B) can change its sign at some value J3* of the mean 



■131 



where the trace r) pp (B = 0) in this equation yields 
Vo"^ = Vt ~ (l/3) e ^ an d ey = e^ej. Therefore, the 



magnetic field [see Eqs. (C17) and (C21) in Appendix 
C]. Thus the point B = B* can determine a steady state 
configuration of the mean magnetic field for e / 0. 
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V. APPLICATIONS: MEAN-FIELD 
EQUATIONS FOR THE THIN-DISK 
AXISYMMETRIC aO— DYNAMO 

Here we apply the obtained results for the nonlinear 
mean electromotive force to the analysis of the thin-disk 
axisymmetric aO-dynamo. Using Eqs. (|30|)-(|32|) we de- 
rive the mean-field equations for the thin-disk axisym- 
metric afi-dynamo: 



dB_ d_f dB\ ^dA 
dt dz \ B dz) dz ' 



OA 



8 2 A 



V *1T 
az 



uB 



(34) 
(35) 



where r, if and z are cylindrical coordinates, B 
Vx(AeJ), G = -r(dtt/dr), and 



Be,, 



Va(B) 
Vb(B) 
V A (B) 

a(B) 



M, 



M, 



M K 
M K 



(Va-Vb)MB\)' 



(v) 



*a(B)ai' 

Lf (1) . 



<J>(B)aW(B) 



and F' = dF/dz, M v = My' + My \ In the axisymmet- 
ric problem dB/dtp — 0. The thin-disk approximation 
implies that the spatial derivatives of the mean magnetic 
field with respect to z are much larger than the deriva- 
tives with respect to r. It is seen from Eqs. (|30|)-(|3^) 
and (^6|)-(^7|) that the contributions to the turbulent dif- 
fusion coefficients tia(B) and tib(B) are from the tensor 
of turbulent diffusion 77^ (B), the tensor ^^(B) and the 
nonlinear velocity U(B) + V^^B). On the other hand, 
contributions to the effective velocity Va(B) are from the 
tensor of turbulent diffusion 7737(B) and the nonlinear ve- 
locity U(B)+VW(B). The functions tja(B), r) B (B) and 
V A (B) are given by Eqs. (|cTl|)-(|C13D in Appendix C. 

The nonlinear dependencies: (A) of the turbulent mag- 
netic diffusion coefficients rjA,(B)/rj^!' and tib(B)/ti^] 
(B) of the effective velocity V A {B)/(B 2 )'; and (C) of 
the nonlinear dynamo number D(B)/D„ are presented 
in Figures 1-3. Here L>* = a*Gh 3 /r]l, D(B) = 
a^(B)Gh 3 /[ VA (B) VB (B)}, 77, = 77^ + (2/3)e^ , a* is 
the maximum value of the hydrodynamic part of the a 
effect, h is the disc thickness and a^(B) = a^$ a (B). 
For simplicity we consider the case e = 0. 

In order to separate the study of the algebraic and 
dynamic nonlinearities we defined the nonlinear dynamo 
number D(B) using only the hydrodynamic part of the 
a effect. We considered three cases: two types of an 
anisotropic background turbulence (ej, = ±1.3577^"'; 
e]L = 0) without magnetic fluctuations (Fig. 1 and Fig. 

3) and an isotropic (e^ = ejf 1 = 0) background turbu- 
lence with equipartition of hydrodynamic and magnetic 
fluctuations (Fig. 2). The negative degree of anisotropy 

implies that the vertical (along axis z) size of turbu- 
lent elements is less than the horizontal size and positive 



(36) 
(37) 



Vb/i (v) t 



0.06 



0.02- 




3 B/B„ 



V< B2 



0.02 



0.01 




0.2 0.4 0.6 0.8 B/B 




FIG. 1: (A). The nonlinear turbulent magnetic diffusion 
coefficients; (B). The nonlinear effective velocity; (C). The 
nonlinear dynamo number for rfc = 0; = — 1.3577^; 



(v) 

Efj, means that the horizontal size is less than the vertical 
size. 

Figures 1-3 and the equations for t)a(B) and t]b(B) 
show that the toroidal and poloidal magnetic fields have 
different nonlinear turbulent magnetic diffusion coeffi- 
cients. In isotropic background turbulence (Fig. 2) the 
nonlinear effective velocity V A (B) is negative. The lat- 
ter implies that it is diamagnetic velocity, which results 
in the field is pushed out from the regions with stronger 
mean magnetic field. In the anisotropic background tur- 
bulence (Fig. 1 and Fig. 3) the nonlinear effective veloc- 
ity is positive, (i.e., paramagnetic velocity which causes 
the magnetic field tends to be concentrated in the regions 
with stronger field). The sign of affects the value of 
r) A (B), t)b(B) and V A (B), e.g., for positive parameter of 
anisotropy the functions t)a(B), ^b(B) and V A (B) are 



= 0. 
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1a.b'i m - 




1 2 3 4 B/B„ 




5 10 15 



B/B. 




V A /(B 2 )' 



DID, 



4 B/B s 




(C) 




B/B„ 



FIG. 2: (A). The nonlinear turbulent magnetic diffusion 
coefficients; (B). The nonlinear effective velocity; (C). The 



nonlinear dynamo number for 77^ 



(v) 



(h) 
Vt ; 



-(") 



= e L h) = o. 



larger at least in one order of magnitude then those for 
negative eft . 

The dependencies of the nonlinear dynamo number 
D(B)/D* on the mean magnetic field B/B cq demon- 
strate that the algebraic nonlinearity alone (i.e., quench- 
ing of both, the nonlinear a effect and the nonlinear tur- 
bulent diffusion coefficients) cannot saturate the growth 
of the mean magnetic field (where B cq — y/JiQpuo). 
Indeed, for anisotropic background turbulence without 
magnetic fluctuations (Fig. 1 and Fig. 3) the nonlin- 
ear dynamo number D(B)/D< t is a nonzero constant for 
(3 ^> 1, i.e., it is independent on (3. This is because for 
(3 ^> 1 the funct io ns t\a oc 1/(3, tjb oc 1/(3 and a oc I//? 2 
[see Eqs. ( |CIS| )-( C2I ) in Appendix C]. In the case of 
isotropic background turbulence with cquipartition of hy- 
drodynamic and magnetic fluctuations (Fig. 2) the non- 
linear dynamo number D(B)/D* oc (3 for > 1 be- 
cause in this case the functions t\a oc 1//3 2 , rjs oc 1/(3 



FIG. 3: (A). The nonlinear turbulent magnetic diffusion 
coefficients; (B). The nonlinear effective velocity; (C). The 
nonlinear dynamo number for rfc = 0; eft = 1.35r)ft ; eft = 
0. 



and a oc I//? 2 [see Eqs. (|Cl^)-(|C2l|) in Appendix C]. 
Note that the saturation of the growth of the mean mag- 
netic field can be achieved when the derivative of the 
nonlinear dynamo number dD(B) / dB < 0. Thus, the al- 
gebraic nonlinearity alone cannot saturate the growth of 
the mean magnetic field. We will show below that the 
combined effect of the algebraic and dynamic nonlincar- 
ities can limit the growth of the mean magnetic field. 
Equation (E8f) in nondimensional form is given by 



da 



(h) 



= 4 



[ VB B'A' - ( VA A" - V A A' 



h 

aB)B] + [C\aft(z)\f v (z)<l> a (B) VA (B)B 2 }' , (38) 



where C is a coefficient, f v (z) describes the inhomogene- 
ity of the turbulent magnetic diffusion, and we define 
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f(z) = a^\z)f n (z). We use here the standard dimen- 
sionless form of the galactic dynamo equation (see, e.g., 
fL6f), in particular, the length is measured in units of the 

disc thickness h, the time is measured in units of h 2 /rj^ 
and B is measured in units of the equipartition energy 
B cq = ytfiopuo- Here uq is the characteristic turbulent 
velocity in the maximum scale Iq of turbulent motions, 
Vt = and and a are measured in units 

of a* (the maximum value of the hydrodynamic part of 
the a effect). For galaxies h/lo ~ 5 and C ~ 0.05 — 0.1. 
Nondimensional equations for A and B are given by 



dB/dt = ( VB B'y - D A' , 
dA/dt = r] A A" - V A A' + aB 



(39) 
(40) 



where Do = a^Gh 3 /^ and B r = —A'{z). In a steady 
state Eqs. ©-© yield 



[MB)B'} 2 + 2CD ^ a (B)7 lA (B)B 2 \f(z)\ = , (41) 

where we used the following boundary conditions B(z = 
±1) = 0, B'(z = 0) = and f{z = 0) = 0. The solution 
of Eq. (Ell) for negative D is given by 



x(B)dB = ^/2C\D^\ VWjdz, 

J\z\ 



(42) 



where x(B) = rj B (B)/[^ a (B)7 lA (B)B 2 } 1 ^ 2 . Consider the 
case e = 0. For f3 > 1 (i.e., for B > l/\/8) the equilib- 
rium mean toroidal magnetic field B(z) is given by 



(43) 



where we used that for )3 > 1 the functions rj A (B) ~ 
3/5/3, 775(B) ~ 2/5/3, $ Q (B) ~ 2//3 2 , and X (B) ~ 2/V£. 



Here for simplicity we considered the case e 



0. 



In a steady state A(z) = — T)b(B)B' /\Dq\. Now we 
specify the profile of the function f(z), e.g., f(z) — 
/4sin(7rz/2)] 2fc+1 [cos(7rz/2)] 2 , where k= 1;2;3;... and 



/* 



2/,' - :-! \ {2k + 3\ {2k+1)/2 
2k + 1 J 



The function f{z) changes in the interval < f(z) < 1 
and it has a maximum f(z — z m ) — 1 at z m = 



(2/tt) arctan[y / (2fc + l)/2]. Equation (|3J) for this profile 
f(z) with k = 2 yields 



B(z) » 0.4C|D |{1 - [sin(^z/2)] 7 / 2 } 2 



(44) 



Equation ( f44|) describes the equilibrium configuration of 
the mean toroidal magnetic field. Thus, the saturation 
of the growth of the mean magnetic field is caused by 
both, the algebraic and dynamic nonlinearities. The dy- 
namic nonlinearity is determined by the dynamic equa- 
tion ( p8| ) whereas the algebraic nonlinearity implies the 
nonlinear dependencies of the turbulent magnetic diffu- 
sion coefficients rj A (B) and t)b{B) and of the effecti ve ve - 
l ocity V A (B) on the mean magnetic field [see Eqs. ( |Cll[ )- 



VI. DISCUSSION 

In this study we calculated the nonlinear tensor of 
turbulent magnetic diffusion, the nonlinear K-tensor, the 
nonlinear effective velocity, and other coefficients defin- 
ing the mean electromotive force for an anisotropic tur- 
bulence. The obtained results were specified for an 
anisotropic background turbulence with one preferential 
direction. We found that the turbulent magnetic diffu- 
sion coefficients for the toroidal and poloidal magnetic 
fields are different. We demonstrated that even for a 
homogeneous turbulence there is the nonlinear effective 
velocity which can be a diamagnetic or paramagnetic ve- 
locity depending on anisotropy of turbulence and level of 
magnetic fluctuations in the background turbulence. The 
diamagnetic velocity implies that the field is pushed out 
from the regions with stronger mean magnetic field, while 
the paramagnetic velocity causes the magnetic field tends 
to be concentrated in the regions with stronger field. 

Note that dependencies of the a-effect, the turbulent 
magnetic diffusion coefficient and the effective drift ve- 
locity on the mean magnetic field for an isotropic tur- 
bulence have been found in |36|, |37], Q using a modified 
second order correlation approximation. Our results are 
different from that obtained in J|6], [ItJ 38 . The rea- 
son is that in |30[ |38| a phenomcnological procedure 
was used. In particular, in the first step of the calcula- 
tions the nonlinear terms in the magnetohydrodynamic 
equations were dropped (which is valid for small hydro- 
dynamic and magnetic Reynolds numbers or in a highly 
conductivity limit and small Strouhle numbers). In the 
next step of the calculations in (3^, [37| Q it was assumed 
that v = 77 = 1 2 /t c , where l c and r c are the correlation 
length and time of turbulent velocity field. The latter 
is valid when the hydrodynamic and magnetic Reynolds 
numbers are of the order of unit. In the present paper we 
use a different procedure (the r-approximation) for large 
hydrodynamic and magnetic Reynolds numbers. 

In this study we also demonstrated an important role 
of two types of nonlinearities (algebraic and dynamic) 
in the mean-field dynamo. The algebraic nonlinearity is 
determined by a nonlinear dependence of the mean elec- 
tromotive force on the mean magnetic field. The dynamic 
nonlinearity is determined by a differential equation for 
the magnetic part of the a-effect. This equation is a con- 
sequence of the conservation of the total magnetic he- 
licity (which includes both, the magnetic helicity of the 
mean magnetic field and the magnetic helicity of small- 
scale magnetic fluctuations). We found that at least for 
the a£l axisymmetric dynamo the algebraic nonlinearity 
alone [i.e., the nonlinear functions a(B), r] A (B), t]b(B) 
and V A (B)] cannot saturate the dynamo generated mean 
magnetic field. The important parameter which char- 
acterizes the algebraic nonlinearity is the nonlinear dy- 
namo number D(B). The saturation of the growth of 
the dynamo generated mean magnetic field by the alge- 
braic nonlinearity alone is possible when the derivative 
dD(B)/dB < 0. We found that for the aSl axisymmetric 
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dynamo the nonlinear dynamo number D(B) is either a 
constant or D(B) oc B for B > B cq /3 depending on the 
model of the background turbulence. Therefore, in this 
case the algebraic nonlincarity alone cannot saturate the 
dynamo generated mean magnetic field. 

The situation is changed when the dynamic nonlincar- 
ity is taken into account. The crucial point is that the 
dynamic equation for the magnetic part of the a-effect 
(i.e., the dynamic nonlinearity) includes the flux of the 
magnetic helicity. Without the flux, the total magnetic 
helicity is conserved locally and the level of the saturated 
mean magnetic field is very low || . The flux of the mag- 
netic helicity results in that the total magnetic helicity 
is not conserved locally because the magnetic helicity of 
small-scale magnetic fluctuations is redistributed by a he- 
licity flux. In this case an integral of the total magnetic 
helicity over the disc is conserved. The equilibrium state 
is given by a balance between magnetic helicity produc- 
tion and magnetic helicity transport Q . These two types 
of the nonlinearities (algebraic and dynamic) results in 
the equilibrium strength of the mean magnetic field is of 
order that of the equipartition field B oq (see Section V) 
in agreement with observations of the galactic magnetic 
fields (see, e.g., Jl6[). 
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APPENDIX A: CALCULATION OF THE MEAN 
ELECTROMOTIVE FORCE 

Let us derive equations for the second moments. For 
this purpose we rewrite Eqs. ([5|) and (Q) in a Fourier 
space and repeat twice the vector multiplication of Eq. 
(||) by the wave vector k . The result is given by 

= ±j2P ip (k)-S ip )SP(b-,B)+S?\b-,B)] 
at flop 

-ft - vk 2 Ul - Fi , (Al) 
= S^(u;B)-S ( t c) (u-B)+G t -r 1 k%,(A2) 

where §!> c) (a;A) = i f a p (k ~ Q)Q p A l (Q)dQ , 
sl b) (a;A) = ik p Ja t (k - Q)A p (Q)dQ , f = 
kx(kxT)/fc 2 , F(k,R,t) = kx(kxF(k,R))/fc 2 p , 
Pij(k) = Sij — kij, 5ij is the Kronecker tensor and 
kij = kikj/k 2 . We use the two-scale approach, i.e., a 



correlation function 

(ui(x)Uj(y)) = J {u i (k 1 )u j (k 2 )) exp{i(k r x 

+k 2 -y)} dki dk 2 
= J fij (k, K) exp (ik-r + iK-R) dk dK , 

= J fij (k, R) exp (ik-r) dk , 

/„(k,R) = J ( Ui (k + K/2K-(-k + K/2)) 
= x exp (iK-R) dK , 

where R = (x + y)/2, r = x y, K = ki+k2, k = 
(ki — k 2 )/2, R and K correspond to the large scales, and 
r and k to the small ones (see, e.g., [B9|, |o|). The others 
second moments have the same form, e.g., 

Mk,R) = y"(Mk + K/2)fo 3 (-k + K/2)) 
x exp (iK-R) dK/pop , 

g tJ (k,R) = /"<&,(]£ + K/2) Ui (-k + K/2)) 
x exp (iK-R) dK . 

The two-scale approach is valid when (1/ B)(dB / dR) <C 
Iq 1 , where B = |B|. Now we derive the equations for 
the correlation functions fij (k, R) , and hij (k, R) , and 
gij(k, R) 



d fy /at 

dhij/dt 
dgij/dt 

i-i i 



i(k-B)$ y + Mij + Ft. 



2vk 2 f l0 , (A3) 
-i(kB)$y + Rij - 2r]k 2 h,j , (A4) 
hj + Cij - {v + rj)k 2 gij , (A5) 



t(k-B)(f K 



hi 



(l/2)(B-V)(^+/ii. 



-fpjBi, p + hi p (2Pji(k) - Sji)Bi iP 

— B p ^qk p (fijq + flijq) , (-^6) 

where V = d/dR, f ijq = (l/2)dfij/dk q , h ljq = 
(l/2)dhij/dk q , and F y (k,R) = (F^k, R) Uj (-k, R)) + 
(w 4 (k,R)F)(-k,R)) , Bi,j = dBi/dRj, and 

(k, R) = [ff tf (k, R) - gji (-k, R)]//i p • (A7) 

The third moments are given by My (k, R) = 
-(f,(k) Uj (-k)) - {ui{k)fj(-k)), R v (k,R) 
{Gi{k)bj{-k)) + (bi{k)Gj(-k)) and Cy(k,R) = 
(^(kK-(-k)) - {kWfji-k)) 



For the derivation of Eqs. (A3)-(A6) we performed sev- 
eral calculations that are similar to the following, which 
arose in computing dgij/dt. The other calculations fol- 
low similar lines and are not given here. Let us define 
Yij (k,R) by 



Y tJ (k,R) = J {sf\u-B-k + K/2)u 3 {-k + K/2)) 

x exp (iK-R) dK = i j \k p + K p /2)B p (Q) exp(iK-R) 
x(tti(k + K/2- Q) Uj {-k + K/2))dKdQ . 



10 



Next, we introduce new variables: ki = k + K/2 Q, 
k 2 = -k + K/2 and k = (k x - k 2 )/2 = k - Q/2, K = 
ki + k 2 = K Q. Therefore, 

Fy(k,R) = i J f t3 (k-Q/2,K-Q)(k p + K p /2)B p (Q) 
x exp (iK-R) dK dQ . (A8) 
Since |Q| <C |k| we use the Taylor expansion 
/y(k - Q/2, K Q) ~ f tJ (k, K Q) 



where ^ = 2(k-B-r) 2 /MoP, % = hkj/k 2 , fl° N) and /^ 0S) 

describe the nonhelical and helical tensors of the back- 

( s) 

ground turbulence. The tensor h\- is determined by an 
evolutionary equation (see, e.g., |3, @, P~5l P3|, p24|, P5|, P^l 



and Section III-C). Now we calculate $^ an d . The 
definition of , given by Eq. ©, and Eq. (|A16|) yield 



l d/ tJ (k,K-Q) 

2 <% s 

and the following identities: 



,y(k,R) pb r(/xop) _1 [^(k.R) - ^(- k : R )] • ( Al ?) 

Substituting Eq. @) into Eq. (|A17|) and using Eqs. (|) 
and (H) we obtain 



Q s + 0(Q 2 ), (A9) 



« {2*(k-B)(/. iJ 



2Bi tP (h pj kn - h lp kij)} .(A18) 



[/y(k, R)B p (R)] K = J fij(k, K - Q)Bp(Q) dQ , Now using Eqs. (^T|)-(|aT5|) and (^T|) we get 
_ r „ /-.. » ~ — ' — ~t~: {2i(k.B)(//^-^f>) 



Vp[/ ii (k,R)fl p (R)] = / i/f p [/ ij (k,R)B p (R)]K 



x exp (iK-R) dK . 
Therefore, Eqs. flA§|)-(|A10|) yield 



(A10) 



F y (k,R) ~ [i(k.B) + (l/2)(B.V)]/«(k,R) 

-k p f ijs (\sL,R)B p>s . (All) 



In Eqs. (A3) and (A4) we neglected the terms oc 



(B-V)gy and oc Bi tP g p j because they contribute to the 
modification of the mean Lorentz force by the turbulence 
effect (see, e.g., |3(], ^l|). In Eq. (A5) we neglected 
the second and hig her derivatives over R. We also ne- 
glected in Eq. (A5) the terms which are of the order of 
Rm _1 V(i?i, fij, hij) and Re~ 1 V(B i ; fy; hij). When the 
mean magnetic field is zero Eq. (|A3|) reads 



M, 



(0) 



R 



(0) 



2vk 2 jf 



(A12) 



We assume that is not changed during the genera- 
tion of the mean magnetic field, i.e., F^ — . This 
implies an assumption of a constant power of the source 
of turbulence. 

Now we split all correlation functions {i.e., 
/,,.//,,.//,,.«l>, ; : into two parts, e.g., f tj = f\P + fjf> 



(N) 



'V 1 '»3 ' 

= [/y(k,R) + / y (-k,R)]/2 and 

/ijV «.,R) — fij(— k, R)]/2. Next, we use r 
approximation which is determined by Eqs. (|g)-( 10 ) . 



where f t 
if } = f.jk.K) 



We assume that r\k 2 <C t _1 and vk 2 <C r -1 for the iner- 
tial range of turbulent fluid flow. We also assume that 
the characteristic time of variation of the mean magnetic 
field B is substantially longer than the correl ation time 
r(fc) for all turbulence scales. Thus, Eqs. ( A3)-([A5|) 
yield 



AN) 
ij 

,w 

AS) 

lij 



AON) 
Jij 

h (0N) 

AOS) 
Jij - 

Tlij , 



(S) 
ij 



-ir(k-B)$ 
- ir(k.B)#|f 
ir(k-B)$W 



(A13) 

(A14) 

(A15) 
(A16) 



(OS) 
ij 



hi?) 



+B (f {0N) 



h 



(ON) 



pi 



) - B hp {f pj 



(ON) 



A N ) 



4> 



(S) 



2zr(k-B) (0 jv) 

[Jij 



(1 + 2V0moP 
+0(B it j) . 



■ (ON), 
l pj > 

(A19) 



(A20) 



To calculate the electromotive force we do not take into 
account the second and higher orders spatial derivatives 
of the mean magnetic field. This implies that in the ten- 
sor bijk (and, therefore, in the tensor $y ) we neglect 
the first and higher orders spati al de rivatives of the m ean 
magnetic field [see Eqs. @, ( |A20[ ) and Eq. flA22| ) be- 
low] . 

Note that Eqs. ( |A13| ) and ( |A14D yield 



(N) 



AN) 



f, 



(ON) 



(ON) 



(A21) 



This is in agreement with that a uniform mean magnetic 
field performs no work on the turbulence. It can only 
redistribute the energy between hydrodynamic fluctua- 
tions and magnetic fluctuation. Analysis in [j3l[ showed 
that a change of the total energy (kinetic and magnetic) 
of fluctuations caused by a nonuniform mean magnetic 
field is of the order ~ rry^ ■* AB 2 . In the case of a very 
strong mean magnetic field Eq. (A21) can be violated. 

Using Eqs. (|A13|)-(|A1^) and (gM-(pg) we calculate 
the electromotive force £ j (r = 0) = J'£j(k) dk, where 

the Fourier component £i(k) — (/iop/2)£imn^lm (k), and 
Eijk is the Levi-Civita tensor. The electromotive force is 



given by Eqs. (|rJ-([L3|). Substituting Eq. (|A14|) into Eq. 
|f) we get 



AON)-, 



-2e imn k mj [h^ k N) - ir(k-B)$^]} dk . (A22) 



Ash 



The integration in k-space in Eq. ( |A22| ) yields 

bijk = £ ijn^nJ (P) + 2£inmCnkjm(P) • 



(A23) 
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Hereafter we use the following definitions: 



X 



(v) 
ijk.. 



X 'l 'J K . . 



x£l \p)-x {h) 



ijk. 



ijk. 

CO 

ijk. 



.(/?) 



and 



A?(/3) = 
C (a) (#) = 



c ij (k)r(fc) 
l + V(/?,k) 
c^(k)r(fc) 
1 + 



dk, 



(A24) 
(A25) 
(A26) 

(A27) 
(A28) 



and Pi = 4B z /{u V2Jw), W,k) = [(/3-k)u r/2] 2 , 

c 



: for Ag> and,, 



/i (0JV) for A W 

For the calculation of the tensor 6yfe we specified a 
model of the background turbulence (i.e., turbulence 
with zero mean magnetic field). The turbulent veloc- 
ity and magnetic fields of the background turbulence are 
determined by Eq . . To i nteg rate over the angles in 
k-space m Eqs. ( |A27^ and ( psj ) we use the following 
identities: 



kij sin 9 
1 + a cos 2 9 

kijmn sin 9 

1 + a cos 2 9 



d9d(p = A 1 5 lj +A 2 p, 



(A29) 



d9 dip — Ci(Sij5 mn + Si m Sj n + Si n Sj m ) 



run mn + ^irnPjn + &inPjm + ^jmPin 

+SjnPirn + 5 m nPij) , (A30) 

where a = [/3u fcr(fc)/2] 2 , ft = ft/ft /3y = /%*>•■ = 
/3;/3,/3 fc . . . , and 



t 2tt 
Ai = — 

a 



A 2 = 



2tt 

a 



arctan (yg) 

(a + 1) = 1 

v a 

arctan(y^) 
(a + 3) -j= 3 



Cl ~ 2^ 



(a+1) 



C 2 = 



Ca = 



2a 2 



2 arctan(i/a) 5a 

arctan(y / a) 



y- 1 



(3a 2 + 30a + 35)- 



7T 

2^ 



(a 2 + 6a + 5) 



arctan(-y/a) 



55a 
13a 



35 



- 5 



To integrate over k in Eqs. ( A27 ) and ( A2S ) we 
use the Kolmogorov spectrum of the background tur- 



bulence, i.e., rfpp N \]i) = r/^'tp^), rhp^'Ck) 



AON), 



ri^ip(k) and ^-° j (k) = (R)(p(k)/3, where ip(k) = 
(nPko^ik/ko)- 7 / 3 , r{k) = 2 To (k/k )- 2 ^ , where k < 
k < kd, ko = Iq , Iq is the maximum scale of turbulent 
motions and kd = fcoRe 3 ^ 4 is determined by the Kol- 
mogorov's viscous scale of turbulence. The integration 



in k-space in Eqs. (A27) and (A28) yields 



^{P) = A^(P)+Pih^(P)P j 



C (a) (P) 



+* 2 (p)n ( fjp n ] , 

+T^\p)5 l3 p m ] , 



(A31) 



(A32) 



where 



^(P)^+^2(P)^P nj 



■hi {[MP) 



l [a) {p) 

r^(p) 

A.a) 



\mp)Wt ] 
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r;;;;,, n = l({[A 2 (p) - cwW P™ 



(A33) 
(A34) 
(A35) 



G 



-C3(/3)(Mt ) /3 J +^/3 P ^ 
^ mn {[A 1 (/3)^ ) + ^C 3 (/3)^ a) ]^ 

+^{^.[C 1 (/3) M W + C 3 (/3)/i^/3 pm ] 
-C^l^Sjn - C 3 (p)fx^>p pm S jn } ,(A37) 



(A36) 



and 



(a) 



= (5/6)[AiC9) + 
A 2 G9) + Ci (/?)], * 2 (/3) = (5/6)[C 3 (/3) - A 2 (/3)], 
* 3 (/3) = [5/3)[A 2 (p)+C 3 {p)]. The functions A„(/3) = 
A n (a)ip(k)k 2 dk = (3/3 4 /tt) ^°(A„(X 2 )/X 5 ) and 
similarly for C n (p), where a = [Pu kT(k)/2] 2 = X 2 = 
P 2 (k/ko) 2 / 3 , and we took into account that the inertial 
range of the turbulence exists in the scales: Id < r < Iq- 
Here the maximum scale of the turbulence Iq <ti Lb, and 
l d = l Q /Re 3/4 is the viscous scale of turbulence, and L b 
is the characteristic scale of variations of the nonuniform 
mean magnetic field. For very large Reynolds numbers 
kd = 17 is very large and the turbulent hydrodynamic 
and magnetic energies are very small in the viscous dis- 
sipative range of the turbulence < r < Id- Thus we 
integrated in A n over k from ko = Iq 1 to oo. The func- 
tions A n (P) and C n (P) are given in Appendix B. In Eqs. 
( A32j )-( [A37 ) we omitted terms which are symmetric in 
indexes i and n because after multiplication Cijm n (P) by 
Sun these symmetric terms vanish [see Eq. (A23)]. 
In order to extract terms cx e 



ijmPm which contribute to 



the nonlinear diamagnetic and paramagnetic velocities, 
we split bijk into two parts, i.e., b 



ijk 



b ijk 



b ifk> where 



b ijk 



P m {S J n[l (P) (P)Pk+^2(P)^ k 1 p 'P P 



-2[T^\p)8 nk p 3 +U^l{P)]} 



(A38) 



12 



*@ = e tjn A^ k >(P) + 2s inm e^ k > m W 



(C) 



(A39) 



[see the definitions given by Eqs. (|A2l)-(|A26|) 1. Next, 
we calculate hj k Bj ik . Using Eqs. (|ll|), ( A38| ) and ( A39| ) 
we also split the electromotive force into two parts 



£ 


= ew+e™ , 


(A40) 


(i) 

i 




(A41) 


(2) 

i 




(A42) 



Using Eqs. ( |A3| ) and ( |A4l| ) we obtain 

ff) = (vWxB) i - 7? g)(VxB), ) 
where 

U i (Ar) (B) = ^[ 7 ^(/3) + 2r( c )(/3)]V l B 2 



(A43) 



1 



+ Jj{ 2U $W + ^W^'PpS^VkB, , (A44) 



B 



(A45) 



and Pij{(3) = — /%. For the calculation of 

the terms oc in these equations we used 

an identity e lmn (i np B m ^ = -[Bx(B-V)B] 2 /B 2 = 
-[BxV(B 2 /2)]i/B 2 - P ip (/3)(VxB) p , which follows 
from the formula (B-V)B = (1/2) VB 2 - Bx(VxB). 
Following to |2l[] we use an identity B±h = (dB)j k — 
£jfc;(VxB);/2 in order to rewrite Eq. ( |A42| ) in the form 



= a.^ + tUxBj.-^tVxB), 

-Kij k (dB)j k , 



(A46) 



where 



(2) 

w 

Kyfc(B) 



(^ P bfl + e jkp bf k > p )/4, (A47) 



(2)' 



(A48) 



-(&gl+©/ 2 



Using Eqs. ( |A43| ) and (|A46 ) we obtain the equation 
for the electromotive force £ = +£^ which is given 
by Eq. (E6L) . The tensor of turbulent magnetic diffusion, 



(A49) 



is given by 

% (B) = S«{AiG9)4 P) 



12 



[C a (/9)+2C 3 (/3)]/# 



+* 2 09) A |f)] + |p 1 +C7 1 )/ lji ] (o) 
+^[C 3 ^] (c) + ^ y [^ 2 (/3)^ P) 



:C 2 (/3)^ 



(^)i 



(A50) 



wher e jtzjj = fi^ (3 nj + f3 in fi $ , a n d we used Eqs. ( |A45| ), 
( |A47| ) and the definitions ( |A24| )-( |A26| ). In particular, 

{X]( c \[3) = JW(|3) - X('")(y2/3) + jrW(V5j8) which 
implies, e. 5 . ; [A 1?7t ]( c ) = A^/3)^ - A x (^)^ + 

Using Eqs. ( A3S|) and ( A4S| ) we calculate Kyfe(B) : 



Kyfe (B) = -5[*i(/3)ig + Mm%\ 

+ b -[{A l - 3d)4- fc ]( c ) - -[OjiV^]^ , (A51) 



r( p )l 



where = s ljn ^ k +e ikn ^J , N$ = Mnp feinfe + 
e ifcn /? pj ) and [C7 3 JV] (C) = C 3 (/3)N$ -C 3 (V20)(N$?j> - 
N^) and similarly for [(Ax - 3Cx)L ijk ]^ [see Eq. 



» 



(|A24|)]. 

The asymptotic formulas for the nonlinear coefficients 
defining the mean electromotive force for /3< 1 are given 
by 



%(B) 



I (M) 
2 M„ 



2 , 
5^ 



0) „W 



42 



(5/i. 



(ft) 



(A52) 



<# } (B) 



ao U) % + - (2/5)/5 2 {^[3a[ ) ' ;) + ^(1 
+(8/7)e)] + ^(2-(9/7)e)}, (A53) 
4 fc) (B)(l-3/3 2 /5)5 y , (A54) 



and for /3 3> 1 they are given by 
Vij (B) 



3tt 
5/3 



(M) 



+^ ) (V2 + 1)])+-|^(9-5V2 
+^ ) (5V2-l)-^[^ ) (5-V2) 



10 1 



+ A W(3 + V2)])-iA 



(-P) 



«& } (B) 



4 ft) (B) 



5 (p) 

■ — [Sij (l-ejvp- (l + 9e)\ 



2a 



(v) 



2/? : 



^2< J (B)% ■ 



(A55) 

(A56) 
(A57) 



The asymptotic formulas for the tensor Kij k for /3 -C 1 
and /? S> 1 are given by Eqs. (^l|) and (I2T 
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APPENDIX B: THE FUNCTIONS A a (/3) AND 
C a (J3) 



The functions A a (j3) and C a (j3) are given by 
arctan (3 



Ai([3) 
A 2 {P) 

C 2 (P) 

cm 



6 
5 

6 

~5 
3 

To 



arctan [3 

p 

arctan /3 



1 



P 



235 
189/3 2 



5 

w 

arctan (3 ( 3 
^ V5 
565 35 
9^ 



5 

15 
+ 7^ 

10 



30 

7/32 



1 

~4 

2 

~ 7 

5 
9/? 1 



5 



(.:-! 



L(/3) 



35 
9^4 



16 
315 



L(f3) 



189/3 2 



arctan (3(1 
5 



6 
7^2 



5 
9^4 



315 



L(J3) 



127 

189/3 2 



5 
9/3 1 



where L(/3) = 1 - 2f3 2 + 2/? 4 ln(l + (3- 2 ). For /3 < 1 these 
functions are given by 



Am 

Ci{p) 
C 3 (P) 



1 



:p 2 



A 2 (p) ~ --/3 2 



C 2 (/3) 



32 
105 



/3 4 ln/? 



35 



/* 2 



and for 3> 1 they are given by 



Ai(/3) 
Ci(/3) 
C3C8) 



3tt 2 

5/3 ~ /J 2 

3tt 

20/3 ' 
3tt 
~20/3 



A 2 (p) 
C 2 (P) 



3?r 4 

'5/3 + p 2 ' 



9?r 
20/3 



Since the function A X {P) + A 2 (/3) ~ 0(/3~ 2 ) for /3 > 1 
(it describes an isotropic part of the a-effect) wc took 
into account in the functions A\{(3) and A 2 (P) the terms 
which are of the order of ~ 0(/3 -2 ). Here we also used 
that for P < 1 the function L(/3) - 1 - 2/3 2 - 4/3 4 In/3, 
and for p > 1 the function L(/3) - 2/3/3 2 . 



APPENDIX C: DERIVATION OF THE 
NONLINEAR DEPENDENCIES 77.4(B), Vb(B) AND 

Va(B) 

Now we consider an anisotropic background turbulence 
with one preferential direction, say along unit vector e, 



where e/3 = 0. In this case 

yW = B - 2 [V^VB 2 + V^e(e-V)B 2 

+V {3) (B-V)B] , (CI) 
U = B- 2 [U (1) VB 2 + U {2) e(e-V)B 2 

+C7 (3) (B-V)B] , (C2) 
Ki jk {dB) jk = -B- 2 {W[VB 2 + 2(B-V)B]xB 

-M«ex(e-V)B}i , (C3) 



-1{A 2 (P)^ T P) + ^C 2 (/3)4 P ) 
+^[(C 2 + 2A 2 ) £Al ]( c )-i[^ T ] (c) } 



7'J 

(C) 



where 
l/< 2 

7J(2 

ir 



In order to derive Eq. ( |C3| ) we used the following identi- 
ties: 

L%l(dB) jk = - e W[ex(e.V)B]i, 



g[(A 2 -C 3 )e p 

-(V2p/48)V{V2p) , 
-(V2/3/4)*' 1 (^/3) £ ( t M ) , 

i* 2 (V2/3)4 M ) , 
-^2(P)e ( -P+5[C 3 e^} 



-^ £ ( l Q )[(Vi? 2 + 2(B.V)B)xB] l 



Using Eqs. (C1)-(C2) and ( A5C) we calculate the func- 



tions M v , M e , Mp, M K , My^ and M<> 2) in Eqs. ( 32 



A!(/3)^ P) + — [Ai(/3) + 4A 2 (/3) + d(/3) 

5 



18 



M K = 
M^ = 

M e = 

M« ee 



/T ' 36 

-C 2 (p) - 5C 3 (p)]e^ 

-[A lVT }^ + ^ 2 {V2P) 
b 

i{3A 2 (/3) ?7 ( p) + [^ 2 ( / 3) + 4vI/ 2 ( / 3)] i 
-* 2 (V2/3)4 M )} , 

-\*i{P)ei P) + \[(C l +A l )e,rK 
+ [/(i) +2W+ i(V (3) + ?7 (3) ) 



[(d+A)^]^) 



■ {M) 



C-P) 
A* 



(C4) 
(C5) 

(C6) 
(C7) 



A 2 (/3)^ P) --[C 2 (/3) + 4C 3 (/3) 



-4A 2 (/3)]4 P ) + i-* 2( V2/3)4 M ) 



14 



~{[A2VT}^ 



[(C 2 

(V2/3/48)*(V2/3) , 



2A 2 



■4C 3 ) £p ] (c) } 



(2) 



^(2) + rj(2 



;[(A 2 -C 3 )e M ] 



(C8) 



(C9) 



Now we take into account that ~V^ N \ U and k contribute 
into the tensor rjij . This implies that in order to calculate 
M n , M e , and Mg we perform the change 

rm -> rjy + P« (/3) [U (3) + C/ (3) + 2W] , (CIO) 



where the second term in (CIO) [which is proportional 
to Pij{[3)] describes a contribution V^ w \ U and n into 
the tensor rjij. Using Eqs. ©-©) and @-(|C9|) we 
calculate the functions tia(B), t)b(B) and Va(B): 



fj(B) + (W/9)[(2C 1 -A 1 )e fl 
5 



(C) 



»/a(B) 

775(B) = ^(B) + ^[(8C 1 + C 2 + 10C 3 -4A 1 



(Cll) 



18 1 

-4A 2 ) £ J^) 



Mi 



A 2 )vT] iC) 



_5 
18 

V2/3 



= {^[(4A 2 - C 2 - 10C 3 )e»]M 



24 



*(v^9)}(ln|B|)' 



(C12) 
[v4 2 r? T ] (C) 
(C13) 



where *(x) = 12[^i(x) + (l/2)A 2 (a;)]^ M) + ^(a:)^* , 
ry(B) = [Ax (/J) + (l/2)A 2 (/3)]^ + (4 P) /12)*o(/U, 
* (/3) = (5/3)[4Ai09) + 3A 2 (/3) + 4Ci(/9) - C 3 (/3)] and 
[X]( c ) is defined by Eq. (A24). The asymptotic formulas 



for the functions tja , ?ys , Va and cU- for p<l are given 



by 



VA(B) = ^-^[3^ + ^(14^) 

= ^-^ 2 rf- 8^ + ^(41^) 



F A (S) = f^ [3j? W_ 4r7 W + 5 (3e ( r) 



-4 £ W)](ln|B|)' 



» 



(B) = <%[(4" 3 -(l/3) £a )(l-^ 2 ) 
2 



6 



and for /3 > 1 they are given by 



(C14) 

) 

(C15) 
(C16) 

(C17) 



MB) = ^[l(V2-l) V ^ + (V2- 7 -)e^ 



(C18) 



Vb(B) = -^-{-[(2V2-l)^ ) + (2V2 + l)4 h) ] 
+ l[(22V2-13)eW + (18>/2 



-i3)4 fe) ]} - 



Vk(-B) 



3tt 

4V2/3 



(ln|S|)' 



(C19) 

(^ ) + ^ ) )+^ ) 
(C20) 



(B) = -^{^ Q e - ^[^ - ~ £q (1 - 6)]} (C21 



„(«) 
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